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Abstract. - This paper is devoted to an analysis of the experiment by Nakamura aZ. (Nature 
398, 786 (1999)) on the quantum state control in Josephson junctions devices. By considering 
the relevant processes involved in the detection of the charge state of the box and a realistic 
description of the gate pulse we are able to analyze some aspects of the experiment (like the 
amplitude of the measurement current) in a quantitative way. 



The possibility to form coherent superpositions of states is one of the most fundamental 
properties that distinguishes quantum from classical physics. On the microscopic level, many 
examples come to mind. Whether it is possible, however, to superpose macroscopically distinct 
quantum states is controversial and has been debated since the advent of quantum mechanics. 
Whereas at the beginning it was thought that the macroscopic world is in some sense classical, 
there have been a number of suggestions over the years on how it might be possible to observe 
macroscopic quantum coherence in solid-state devices Q]. Superconducting nanocircuits have 
been used successfully to test quantum mechanics in mesoscopic systems. Examples are the 
test of the Heisenberg uncertainty principle in a mesoscopic superconductor Q| or experiments 
on the superposition of charge states in Josephson junctions Q . More recently, the increasing 
interest in quantum computation |Q and the search for implementations that can be scaled 
and integrated has made superconducting circuits promising candidates to realize qubits 

In a recent experimental breakthrough, Nakamura et al. Q demonstrated coherent oscil- 
lations between two charge states of a superconducting island in a single-electron device, the 
so-called Cooper-pair box. These states are macroscopically distinct in the sense that they 
are two different states of an island that contains, say, 10* electrons. The superposition is 
achieved by switching the gate voltage of the box quasi-instantaneously to the point at which 
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Fig. 1. - Equivalent circuit. The Cooper-pair box corresponds to the configuration Jl = E.j, Jr = 0, 
cVl = and cVr = —eV ~ 2A using the right lead as a probe gate. 



the two charge states, n — and n ~ 2, are degenerate. Here n is the number of excess charges 
on the island (one excess Cooper pair means n — 2). As a result, the system performs Rabi 
oscillations between these two states that are monitored by another weakly coupled tunnel 
junction. 

In this paper, we analyze the experiment by Nakamura et al. Q by solving the appropriate 
master equation. The reason for performing the work presented here is twofold. Although the 
description of the coherent oscillations in terms of two-state system is appealing and contains 
most of the physics there are fundamental questions which are still unanswered. They are 
related to the mechanisms of decoherence and, particularly important for this experiment, the 
role of the measuring apparatus and real shape of the gate pulse. Each of these issues leads 
to computational errors when the Cooper-pair box is used as a (charge)-qubit. Therefore to 
investigate each of these issues in great detail is a necessary prerequisite for the implementation 
of a solid state quantum computer with Josephson nanocircuits. The next step in the coherent 
control of the dynamics of a macroscopic system is the experimental verification of conditional 
dynamics, i.e. of entanglement. Entangled states, a central concept of quantum mechanics, 
are difficult to characterize and measure. A quantitative understanding of the single qubit 
experiment is also relevant for modeling two-qubit gates. 

To investigate the results obtained in Ref. |Q, we consider the superconducting transistor 
shown in Fig. |l| and described by the Hamiltonian (see, e.g., [^) 

The first two terms define the charging part (He), for which we have adopted the effective 
capacitance model: C = Cl + Cr + Cg with junction capacitance Cl and Cr for the left and 
right junction and Cg for the gate. When the leads and the gate are biased by voltages Vl, Vr, 
and V^, respectively, the total gate-induced charge is given by Qt = CLVL+CRVR+CgVg. Qj/e 
is the number of electrons which have passed through the junction j to the central electrode 
D, and Q — Ql + Qr. is the total charge on D. The last term in ([^) is the Josephson part 
(iJj) with the Josephson coupling energy Jj and the phase difference (pj across the junctions 
j = L,R. In the Coulomb-blockade regime, where the charging energy Ec — e^/2C is larger 
than Jj , the charge Q on the island is quantized in units of electric charge e. Accordingly, we 
will use the basis of charge states \n) with charge Q = ne where n is integer. 

The coupling of the system described by Eq. to the environment leads to decoherence. 
The environment or bath that we focus on here are the quasiparticles on the two supercon- 
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ducting leads L, R and the central island D 



(2) 



a=L.R,D ka 



Here, 7]^^ ilka) creates (annihilates) a quasiparticle with momentum k and energy e^a — 
\/ if.^ + in electrode a. Also, i^fc is the usual single-particle dispersion (with respect to 
the chemical potential), A is the superconducting gap (to simplify the discussion, we assume 
identical superconductors, as bulk, for all electrodes) and a labels the spin. The effects of 
quasiparticle tunneling can be described by the tunneling Hamiltonian 



(3) 



where Xj = X^fc^o- '^kq^l-j^JgDa- ^^'^ "^kq is the tunneling amplitude. The total Hamiltonian is 
given by Htot ^ H + Hgp + Ht- 

The evolution in time of the reduced density matrix p — Tr^i ,yPtotai is governed by the 
generalized master equation (see, e.g., §) 



h 

E 



[H,p(t)] 



ds af{s) 



E / dsaf{~s) 



Pit) 



h.c. 
- h.c. 



(4) 



The correlation functions 



■{t) EE (^X,{t)x]) Ih^ and a<{t) ee (^x]x,{t)) /h' 



describe the 



effects of quasiparticle tunneling. More explicitly, the Fourier transforms of these correlation 
functions can be expressed in terms of the quasiparticle current-voltage characteristics /J^(£'): 
a>{E) = [l+nB{E)]I]P{E)/e and afiE) = nB(-B)/f (-B)/e, where ngiE) = 1/ (e/^^ - l) is 
the thermal distribution function. 

The contribution to the transport due to resonant Cooper-pair tunneling has been consid- 
ered by Averin and Aleshkin |]lOj and by van den Brink et al. ||l^ . They were interested in the 
d.c. transport current. In the present paper, we focus on the coherent oscillation between the 
charge states and its decoherence due to quasiparticle tunneling. In this case, the off-diagonal 
elements play a crucial role in the dynamics of the reduced density matrix. We note in passing 
that it is also possible to use He instead of H as unperturbed Hamiltonian and the 
Josephson tunneling Hj in (|^) and the quasiparticle tunneling Ht in (^) as perturbations; 
this is equivalent to our model to second order in the tunneling amplitude. 

The Cooper-pair box is probed via quasiparticle tunneling |7j. This can be achieved by 
configuring the system parameters as follows: the whole system is biased so that Cooper-pair 
tunneling occurs only across the left junction whereas there is only quasiparticle tunneling 
across the right junction, cVl — and cVr = —eV ^ — 2A. Respectively, we can safely set 
Jl = Ej and Jr — 0, and it is always implied that ksT <^ Ej <^ Ec <S A. Due to the strong 
Coulomb repulsion Ej <C Ec, it suffices to consider the subspace of {|0), |2)}. In this 
basis, the time evolution of the diagonal elements a-„ = p„^„ (n = 0, 1, 2) and the off-diagonal 
element x = P0.2 is given by 



x1-r+(o)ao(i) + r-(i)ai(0 



(5) 
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(72(0 = +*§[x-x1-r-(2)a2W+r+(l)ai(i) (6) 
= — — x{t)-t^[<Jo-cr2] (7) 

-i[r+(o) + r-(2)] xW, 

together with the normahzation condition Trp[t) = ao{t) + cri(t) + a2{t) = 1. Here £^2,0 — 
4[1 — Qt(i)/e]-Ec is the change in the charging energy when a Cooper pair tunnels into the 
Cooper-pair box across the left junction. In Eqs. (^), r^(n) = rf^{n) + T^{n) and ^fin) is the 
quasiparticle tunneling rate across the junction j = R resulting in the transition n — > n ± 1. 
These rates can be written explicitly as 

r+(0) = a>{E^)+a<{eV-E^) (8) 

r-(2) = a>{E+)+a>{eV + E+) (9) 

r+(l) = a<iE+) + a<ieV + E+) (10) 

r-(l) = a<{E^)+a>{eV-E^) , (11) 

where E± — Eq [1 ± 2(1 — Qt/e)]. It is useful to notice that at sufficiently low temperatures 
(fcsT < A), r+(0) and r+(l) are exponentially small (~ e'^/'^'^^) while r-(2) ~ r-(l). 
From (|[), it is clear that quasiparticle tunneling causes decoherence in the system. Indeed, at 
Qt = e, the coherent oscillation between the two degenerate charge states |0) and |2) decays 
with the time scale l/F = l/r-(2) ~ l/r-(l), 

'^oit) ^ -H \ — icosut + —sin ujtj (12) 

I g-3rt/2 e"'"*/^ 
a2{t) ~ o+- — 2 o — cosujt (13) 
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3-Bj 12^j 



2sintjt cosijjt , (14) 

J 

where the oscillation frequency is also slightly modified as w ~ (Ej/h) \l + (hr/Ejy\. Here 

we assumed the initial conditions cro(O) — 1 and cr2(0) = x(0) ~ 0. In the stationary limit 
{t 00), the charge states 0), |1), and |2) are equally populated. This corresponds to the 
Josephson-quasiparticle cycle |ll|, |l2|, i.e., resonant tunneling of Cooper pairs followed by 
sequential tunneling of quasiparticles to the probe gate. 

We now look more closely at the experiment by Nakamura et aL |Q. A pulse of finite 
length At was applied to the gate to change the total gate-induced charge from Qt = Qo 
to Qt — Qo + Qp and back. For example, if Qo/e < 1 (far from the resonance point) and 
(Qo + Qp)/c- ^ 1 (s^t the resonance point), the pulse causes the system (which is initially in 
state |0)) to oscillate between |0) and |2). Depending on the pulse length At, the system may 
or may not be in |2) with an increased probability at the end of the pulse. The decay of |2) 
to |0) through quasiparticle tunneling leads to an excess current. This was illustrated in 
by solving the time-dependent Schrodinger equation for an isolated two-level system. Here we 
solve the master equation (^, which includes the decoherence due to the measuring device, 
and calculate the time-averaged current / = {I{t)) = r~ {2){a2{t)) +r~ {l){ai{t)) . To simulate 
the experimental situation, we generated an array of pulses with repetition time Tr (after a few 
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Fig. 2. - (a) Plot of pulse-induced current versus Qo /e and At. (b) Current through the probe junction 
(R) versus Qo/e with (solid) and without (dashed line) pulse of length At = 2Tcoh — 160ps. In both 
plots, we used the resistance ratio Rl/Rr = 1/1800 which corresponds to quasiparticle tunneling 
rates r-(2) - (6ns)-i and r-(l) ~ (8ns)~\ i.e., like in 



times of repetition, the time series of p{t) repeats the same pattern). We took the parameters 
used in the experiment, i.e., Ej — 51.8/^eV, Ec — 117/xeV, and A = 230/ieV. The voltage 
was chosen to be eV 2A + 1.65Ec, and the repetition time Tr = 16ns. Figure ^ shows the 
result of our calculation, viz., the pulse-induced excess current versus offset gate charge Qo/e 
and/or pulse length At. 

Apart from decoherence by quasiparticle tunneling, the real experimental situation includes 
several complications: (i) the repetition time is finite and comparable to l/r~(2) and 
l/r^(l); (ii) 'jittering' of the pulse; and (iii) finite rising/falling time of the pulse. The effect 
of these complications is summarized in Fig. ^ and they are discussed below. 

Ideally, the time between two pulses should be long enough (compared to l/r~(2) or 
l/r~(l)) such that the system is guaranteed to relax to the desired initial state. Experimen- 
tally, Tr is limited by the detector sensitivity, since the maximum value of the pulse-induced 
current is given by A/max = 2e/Tr. The value of Tr in was comparable to l/r~(2) or 
l/r^(l) and the waiting time, T^ =Tr — At, was not long enough. First of all, this reduces 
the pulse-induced current collected in the probe gate with respect to the ideal maximum current 
A/max- And secondly, the finite Tr also results in a wiggly behavior of AI{At): The charge 
oscillations with small amplitude and large frequency (~ \/ E'^ + E^j/h) have not decayed yet 
when the next pulse is turned on, see the time dependence of (T2 {t) in Fig- (a) (dashed line) . 
A small change in At with Tr fixed leads to a change in T^,, which may be comparable to 
h/ ^ E'^ + Ej and hence results in a rapid change in the pulse-induced current, see Fig. || (b). 

The jittering of the pulse leads to fluctuations in Tr (or T^,). This smears out the wiggly 
behavior in AI discussed above. The pulse jittering in |^] was of the order of several percent 
of the coherence time Tcoh — h/Ej [|l3j. Numerically, the jittering can be taken into account 
by averaging AJ(Ai) over slightly different values of At. In Fig. || and Fig. || (b) (fat solid 
line), we took the average over the interval (Tr — 7ps, Tr + 7ps). 

Experimentally, the oscillation amplitude of the pulse-induced current measured in the 
probe gate was significantly smaller than the ideal maximum value A/max- This may be 
accounted for by taking into account the finite repetition time and finite rising/falling time of 
the pulse. The effect of a finite repetition time was already discussed above. To understand 
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Fig. 3. - (a) Pulse shapes used to create the current in (b) and corresponding time dependence of (J2(i). 
Each step of the step-wise pulse is 0.2Tcoh long. The other parameters used are Tr = 200Tcoh — 16ns, 
Qo/e = 0.51, Rl/Rr = 1/1800 (corresponding to the values of r"(2) ~ (6ns)~^ and r~(l) ~ (8ns)"^ 
in |7|). (b) Pulse-induced current versus pulse time At with finite (thin solid line) and zero (dashed 
line) rising/falling time of the pulse. The fat solid line shows the effect of additional averaging by 
jittering of the pulse; its amplitude reproduces the experimental value. 



the effect of the finite ramping time of the pulse, we increased/decreased the gate voltage in a 
step- wise way, see Fig. ^ (a). As sliown in Fig. || (b), the oscillation amplitude for this pulse 
shape (thin solid line) is suppressed compared with that for an ideal pulse (dashed line). The 
fat solid line shows the effect of additional averaging by jittering of the pulse; its amplitude 
reproduces the experimental value. The shift of the oscillation in At is due to the change in 
the effective pulse length. 

In conclusion, we have presented an analysis of the experiment by Nakamura et al. 
by solving the appropriate master equation. In particular, we have considered the relevant 
processes involved in the detection of the charge state of the box and have used a realistic 
description of the gate pulse. 

We would like to acknowledge stimulating discussions with Y. Nakamura, G. Falci, and E. 
Paladino. 
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